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We investigate the tunneling process between two sym- 
metric stable islands of a forced pendulum Hamiltonian in the 
weak chaos regime. We show that when the tunneling doublet 
is quantized over a classical non-linear resonance the tunnel- 
ing rate strongly deviates from the semiclassical prediction. 
This mechanism is responsible for the irregular dependence 
of the tunneling rate on the system parameters. The weak- 
chaos condition allows us to make a theoretical prediction 
which agrees very well with the numerical results. This opens 
up a possible avenue to a general theory on the dependence 
of quantum tunneling on classical chaos. 
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One of the most attractive aspects of classically chaotic 
quantum systems is the surprising behavior of the tun- 
neling rate. This is shown to be an erratic function of the 
system parameters This means that the tunneling 

rate undergoes fluctuations which enhance, or reduce, its 
intensity by several orders of magnitude, compared to 
the smooth and regular behavior corresponding to the 
traditional conditions. 

A kind of general agreement has been reached by the 
researchers working on this hot issue. This is that the 
tunneling properties must be traced back to the cross- 
ing of the tunneling doublet with a third state Q| , which 
corresponds to the chaotic region of the classical phase 
space Q. We shall be referring to these phenomena as 
Chaos Assisted Tunneling (CAT) processes, since this is 
the definition generally adopted in literature in spite of 
the fact that the influence of chaos can also reduce the 
rate of the tunnel process. It is also thought, at the level 
of merely qualitative arguments, that the erratic-like be- 
havior of the tunneling rate depends on the chaotic na- 
ture of the third state. It has to be pointed out that a 
semiclassical treatment represents the natural way of es- 
tablishing a connection between the quantum properties, 
tunneling rate in the case under discussion, and the clas- 
sical phase space. Unfortunately, these methods cannot 
be applied to the chaotic states, and this is probably the 
reason why all the results reached so far on this inter- 
esting issue fail supplementing the qualitative arguments 
with precise quantitative predictions. 

This letter is devoted to describing the discovery that 
the processes observed in literature so far, the CAT pro- 



cesses, are a special case of a more general phenomenon. 
The third state invoked by all the researchers of this field 
need not be chaotic. Here we show that the same qualita- 
tive behavior as that observed in [|],|| becomes ostensible 
also as an effect of the crossing with regular states. In 
this case the crossings of states, and thus the fluctua- 
tion of the tunneling doublets, can be directly related to 
the classical phase-space structure and in particular to 
the destruction of the regular tori by the non-linear res- 
onances. This makes it possible for us to derive a quan- 
titative prediction, which turns out to agree remarkably 
well with the numerical results. 
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FIG. 1. Stroboscopic plot of the classical dynamics driven 
by Hamiltonian (|l]) . The values of the system parameters are 
Q = 2, vi = 0.035, i>2 = 0.005. Some non-linear resonances 
of order 1/n are visible: those below the separatrix refer to 
n=6, 7 and 8, and that above it to n~8. 

To address the problem of the connection between tun- 
neling and classical phase-space structure we consider the 
following Hamiltonian 

V 2 

H = y + vi cos(2g) + v 2 cos(2q - fit) . (1) 

This is a forced system whose dynamical properties are 
well known. Depending on the value of the perturbation 
strength v 2 , the phase space can show either regular or 
chaotic dynamics. For v 2 = this is nothing but an or- 
dinary, and integrable, pendulum. Increasing the value 
of v 2 makes the dynamics non integrable and eventually 
chaotic. Fig. [j] shows the effect of setting v 2 = 0.005. 
We can easily identify the isolated resonances of order 
1/8, 1/7 and 1/6, resulting from the destruction of tori 
with winding numbers rationally connected to the pertur- 
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bation frequency. Notice also the stochastic layer close to 
the separatrix, where the higher-order resonances over- 
lap, giving rise to fully developed chaos. 

To study the quantum dynamics of this system, and in 
particular its tunneling properties, we set periodic con- 
ditions at the borders of the interval q £ [0, 2ir]. This 
has the effect of creating an enlarged phase space con- 
sisting of two identical cells. The quantum mechanical 
eigenstates must be either symmetric or anti-symmetric 
with respect to a 7r-translation along the g-axis, and some 
eigenstates exist whose linear combination is essentially 
located in only one of the two cells. According to || 
the time evolution is described by means of a unitary 
operator, which evolves the quantum state by one entire 
period of the external perturbation. This is called Flo- 
quet operator and it is denoted by F. The generalized 
eigenstates u n (q,t) and eigenvalues e„ of F are derived 
from the eigenvalue equation: 

d 

(H - i?i—)u n (q, t) = e n u n (q, t) , (2) 
ot 

where u n (q,t) — exp(— i(/) n t)^0 n (g, i) is by construction 
time periodic, u n (q, t) — u n (q, t+T), and e n is defined by 
£„ = h(/) n . One important feature of the Floquet formal- 
ism is the Brillouin-zone structure of the spectrum: every 
eigenstate u n (q, t) results in the wider class of eigenstates 
and eigenvalues defined respectively by 

u n , m {q,t) = u n (q,ty mnt (3) 

and 

£n,m = £n + ^mVl , (4) 

where rn e Z. It is evident that the functions belonging 
to the same class represent the same solution of Eq. (^) . 
Therefore, if a eigenstate has a given eigenvalue e, which 
does not belong to the first Brillouin region, < e < fi,f2, 
it is convenient to fold it back to this region by means of 
the prescription (Q). To avoid confusion we denote the 
resulting energies with the symbol E. 

The eigenstates u n}Tn (q,t) reflect the translational in- 
variance q — » q + it and can be labelled as odd or even 
with respect to this symmetry operation. Thus the en- 
ergy spectrum of the Floquet operator appears as a series 
of doublets E n} ± with ± indicating the translation sim- 
metry. The rate of the tunneling process is determined 
by the energy splittings AE n = \E n<+ — 2? n ,_| which, in 
turn, determine the dynamics of a generic wavepackct 
initially located in one of the two islands. For simplic- 
ity here we focus on the energy splitting of the "main 
doublet" , which represents the natural extension of the 
fundamental-state doublet of an autonomous Hamilto- 
nian to the Floquet picture. This is defined as the set of 
the two eigenstates of the Floquet operator that have the 
largest overlap with a minimum-uncertainty state located 
at the center of one of the two stable islands || . 

We adopt the following numerical procedure. The Flo- 
quet matrix F is determined by the numerical integration 



of (g) and the eigenvalues and the eigenstates are sub- 
sequently obtained by numerical diagonalization. The 
results for the main-doublet splitting AE = AEq, as a 
function of v\ are shown in Fig. |^a for different values of 
the perturbation parameter v-i- We see that accordingly 
with the semiclassical prediction in the unperturbed 
case AE is a smooth function of v\. However, when 
V2 is given an even small but finite value, AE strongly 
departs from the smooth behavior, with an increase, or 
decrease, by several orders of magnitude. This behavior 
is similar to that of 0, where the splitting irregularity 
has been related to the crossings of the tunneling doublet 
with a third level which belongs to the chaotic region of 
the corresponding classical system. To confirm the con- 
nection of these results with the third-state theory, we 
calculate the first energy levels E n . In Fig. ||b we show 
those of them which cross the main doublet in the same 
Vi— interval as that of Fig. ||a. We see clearly that the 
splitting irregularities correspond to the crossing between 
the main doublet and a third state, or more precisely, a 
second doublet p|. The effect of level crossing with a 
weakly perturbed Hamiltonian is well established: the 
spectrum is modified only in the vicinity of the crossing 
which becomes an avoided level crossing, thereby result- 
ing in significant changes of the main-doublet levels, and 
thus of the splitting AE. Note that the avoided crossing 
is too small to be visible in the scale of Fig fy. 

Is this effect a CAT process? Do the states colliding 
with the main doublet belong to the chaotic region of the 
phase space? The answer is incontrovertible: in general 
they do not. This is so because, at least for the smallest 
values of v-i used in Fig. ||a, the doublets 5, 6 and 7 be- 
long to the regular region. In fact v-i = 0.0001, a value at 
which the fluctuations of the tunneling rate are already 
very strong, means a perturbation weaker than that used 
to derive Fig. [I]: the corresponding phase space is even 
more regular, and thus barely distinguishable from the 
unperturbed one except possibly for a thin stochastic 
layer around the separatrix. We think, therefore, that 
the results of Fig ^a, with the help of Fig. [j], afford a 
compelling numerical evidence that the tunneling rate 
can be characterized by strong fluctuations without nec- 
essarily involving the interaction with a strongly chaotic 
region. The splitting irregularities are shown below to be 
caused by the birth, in the regular phase-space region, of 
non-linear resonances, a fact that, surprisingly enough, 
has been overlooked by the literature on this field of re- 
search. 

Note that further progress on the CAT processes has 
been hampered by a major difficulty: it is not yet known 
how to realize a proper semiclassical picture of the chaotic 
states and, consequently, how to evaluate the tunneling 
matrix elements. The crossings 5, 6 and 7 of Fig. |^a, on 
the contrary, are compatible with the semiclassical quan- 
tization of Hamiltonian ([l]) , along the lines established by 
Breuer and Holthaus 0. The analysis here is restricted 
to the crossing between the main doublet and the states 
lying below the separatrix. 
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We proceed as follows. Adopting a perturbative ap- 
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proach, we replace H with Ho = \ + UiCOs(2g), and 
write the condition of level crossing in the Brillouin zone 



H Q {K{n + i)) + nClm = H {h{n' + i)) + Mm'. (5) 

Let us assume h to be so small as to make it possible to 
consider Ti(n — n') to be a small expansion parameter. By 
using J = h(n + 1/2), we obtain 



Am 
An 



Q + 0(h) , 



(G) 



where U = u>(J) — ^f- is the frequency of the unper- 
turbed libration as a function of the classical action J, 
An = n' — n and Am = m — m'. Eq. (^|) corresponds 
to the classical condition for the onset of non-linear res- 
onances. This means that, for sufficiently small h's, the 
crossing of the two unperturbed levels E® m and E°, m , 
occurs as a quantum reflection of the birth of a non-linear 
resonance of order Am/ An in the classical phase space. 
The latter has the same energy as the torus determining 
the semiclassical quantization of one of the two crossing 
levels. On the basis of this result the level crossing con- 
cerning the main doublet (n = 0) corresponds to the over- 
lap between the semiclassical quantization torus of the 
fundamental state and a non-linear resonance of the ap- 
propriate order. Notice that, in principle, Eq. (^) involves 
resonances of any order thereby making a given level un- 
dergo a virtually infinite number of crossings upon change 
of a system parameter. However, Ref. fij shows that the 
perturbation strength at the crossing is determined by 
the order of the corresponding non-linear resonance, and 
that the first-order resonances (Am = 1) produce the 
most intense effects. 

We have now to point out that our numerical results 
are confined to finite, and relatively large, values of Ti by 
the obvious limitations of computer calculations. This 
does not conflict with physical reality, since an exper- 
imentally observable tunnel process implies the Planck 
constant to be finite rather than infinitelly small. In con- 
clusion, to derive a more realistic prediction we must go 
beyond Eq. (|^), and this, at the same time, makes it pos- 
sible, as we shall see, a reliable comparison between the 
semiclassical prediction and the numerical results. This 
purpose is realized by applying to Eq. (||) the second- 
order approximation. We obtain: 



Am 



n 



An (i + |||An) 



(7) 



where 



4§ can be expressed analitically via uj = 
TT^/vi/K(k), K(k) is the elliptic function and k 2 = 
(E + v 1 )/2vi. Eq. (@) looks like Eq. (§) with the fre- 
quency properly renormalized and thus depending on 
h, as well as on the energy of the colliding level. 




FIG. 2. The origin of the tunneling rate fluctuations for 
Hamiltonian (|l|) with Q — 2 and Ti — 0.025. (a) The main 
doublet splitting AE as a function of vi. The thick full line 
denotes the unperturbed case (v2 = 0). The other curves, 
from the bottom to the top, refer to v% = 0.0001, 0.001 and 
0.005. (b) The levels E„ as functions of v\. The full line la- 
belled by 0± denotes the main doublet. The other levels are 
indicated by either full or dashed lines according to whether 
their crossing with the main doublet occurs below or above the 
separatrix. The full lines are labelled by their quantum num- 
bers. The value of the perturbation strength is «2 = 0.001. 
All the energy splittings, also at the avoided crossings, are 
not visible on this scale, (c) The energy of the fundamen- 
tal doublet (full line) and the energies of the renormalized 
quantum resonances as solutions of Eq. (Q) (symbols). The 
resonances are of order 1/n = 1/5 (circles), 1/6 (squares) and 
1/7 (triangles). The vertical arrows are eye-guides connecting 
the theoretical crossings of (c) to the numerical crossings of 
(b), and these to the tunneling peaks of (a). 

This theoretical prediction on the quantum crossings 
is illustrated by Fig. ^p. We note that the fundamental 
doublet energy is denoted by a single full line because, 
as in Fig. ||b, the energy splitting is not visible in this 
scale. The prediction stemming from Eq. (^) must be 
compared to the numerical quantum crossings of Fig. ||b, 
which occur both above and below the separatrix. Here 
we see that the accuracy of the prediction is good and, 
as expected, becomes better with resonances of smaller 
order. A still better agreement is expected from smaller 
values of h. 

We now invite the reader to compare ||b to ||a, while 
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keeping in mind its earlier correspondence with gc. This 
results in a vivid explanation of the "fluctuations" of the 
tunneling rate. This is made especially evident by a thor- 
ough examination moving from smaller to larger values of 
the perturbation strength V2- We see that at i>2 = 0.0001 
the deviations from the unperturbed behavior are signif- 
icant only in the close vicinity of the level crossings, thus 
lending support to our perturbative approach. With the 
increase of v%, the peaks broaden and overlap one an- 
other and for still larger values of t>2 new peaks come 
forth corresponding to crossings of higher order. 

The peaks around Vx = 0.025 and V\ — 0.0275 refer to 
levels close to the separatrix, which at properly large v's 
are embedded in the chaotic region. The peak at about 
v i = 0.020 refers to an above-separatrix level. As a con- 
sequence of this, the dynamical classical process behind 
these three peaks is compatible with the motion from the 
one to the other potential well. For these reasons one 
would be tempted to conclude that the corresponding 
transition rates are the largest. We see, on the contrary, 
that the peaks referring to a below-separatrix level cross- 
ing (v2 > 0.030) are more intense than the peaks due to 
the above-separatrix level crossings. We also note that 
at the highest value of v 2 used in Fig. |]a, when the tun- 
neling rate results in the maximum departure from the 
unperturbed prediction, the classical phase-space struc- 
ture is still almost regular, as shown by Fig. [I], which 
refers to the same value of «2- 

In conclusion, we have identified a process resting on 
the role of isolated non-linear resonances rather than that 
of the connected chaotic sea. The corresponding tunnel- 
ing rate AE exhibits fluctuations which can be even more 
intense than those provoked by the crossing with chaotic 
states. This can be accounted for by using the same 
heuristic arguments as those adopted in Ref. ||. These 
authors argue that the intensity of the process is pro- 
portional to the matrix element V connecting one of the 
partner of the tunneling doublet, state \r), to the crossing 
state |c). This crucial parameter is roughly proportional 
to the overlap between \c) and |r). We can predict its 
value with heuristic arguments based on the observation 
of the Husimi distribution of the eigenstates (not shown 
here). The ground state is represented by almost gaus- 
sian packets located at the centers of the regular islands, 
the regular crossing state corresponds to a double ring- 
shaped bun surrounding the gaussian packets, and the 
chaotic state to an eight-shaped distribution lying on the 
chaotic region of Fig. [I]. As the value of % decreases, the 
gaussian state shrinks, the regular distribution collapses 
on its quantization torus and the chaotic state remains 
localized in the chaotic layer. This produces a decrease 
of the overlap among different states, and consequently 
an overall decrease of the tunneling rate. However, a 
mere inspection of Fig. [j] leads to the incontrovertible 
conclusion that the effect of the regular-regular crossing 
is expected to remain larger than that of the regular- 
chaotic one. This is made compelling by the topology of 
the phase space: the overlap with the closer regular state 



remains larger than the overlap with the farther chaotic 
one. 

We note also that the values of the parameter v\ cor- 
responding to the tunncling-rate peaks, appear to be in- 
dependent of the intensity of the perturbation strength 
i>2- Our conjecture is that the position of the peaks re- 
mains unchanged even when the overlap of the isolated 
resonances occurs, and a condition of full chaos devel- 
ops. This changed condition produces the broadening 
of these peaks, the birth of higher-order crossings and 
the resulting merging of all of them into a single erratic- 
like structure, without influencing their original position. 
Thus we think that even after the crossover to the chaotic 
regime, significant signs of the peaks originated by the 
quasi-integrable non-linear resonances remain. 

Before concluding, we want to point out that these 
results are qualitatively similar to those obtained in 
Ref. These physical effects, as here shown, can be 
accounted for very well, both qualitatively and quantita- 
tively. This means that the theoretical arguments of this 
letter provides an exhaustive account for the "plateau 
effect" of This is an important task which has re- 
cently challenged the efforts of some groups (see Ref. M ) . 
Should the phenomenon illustrated in Q be proved to be 
a form of CAT, this would be a strong support of our 
conviction that the semiclassical arguments of this let- 
ter are a significant step ahead towards a quantitative 
prediction of CAT processes. 
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